Flow-Equation method for a superconductor with 

magnetic correlations 



C.P.Moca 
Dept. of Physics, Univ. of Oradea 
M.Crisan and I.Tifrea 
Dept. of Theoretical Physics, " Babes-Bolyai" University, Cluj 

February 1, 2008 

Abstract 

The flow equation method has been used to calculate the energy of single impurity 
in a superconductor for the Anderson model with U ^ O.We showed that the energy of 
the impurity depends only of the (renormalized order parameter) which depends of 
the renormalized Hubbard repulsion . For a strong Hubbard repulsion = U and 
A-^ = A^ the effect of the s — d interactions are nonrelevant,a result which is expected 
for this model 
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1 Introduction 



The flow equation method given by Wegner Q has been successfully applied for the 
many-body problem by Kehrein and Mielke 0, for the Anderson Hamiltonian. In 
a previous paper the present authors |Q] showed that this method can be used to 
calculate the energy of a superconductor containing magnetic impurities describe by 
the Anderson Hamiltonian. with U = 0. We studied (See ref.|Q])the influence of the 
density of states on the single impurity energy for the case of a van-Hove density of 
states. In this case the energy is reduced by the superconducting state and corrections 
depends on A^. In this paper we consider a superconductor with a constant density of 
states but for the impurity we take the Hubbard repulsion U ^ 0. 

2 Model 

We consider a superconductor containing magnetic impurities describe by the model 
Hamiltonian. : 

H = Hbcs + Ha (1) 

where Hbcs is 

Hbcs = ^fc^aCfc.^ + ^(^t^ifca + Cfc,TC-fca) (2) 

fc,(T k 

A being the order parameter and Ha the Anderson Hamiltonian. describing the impu- 
rity in a metal, as: 

Ha = Y ^ddtda + VkA4,ada + d+Cfc,,) + Ud+d+d^d^ (3) 
d k,cr 

In equation (3) the first term is the energy of the impurities, the second term is the 
interaction between the itinerant-electron and impurity and the last term in the Hub- 
bard repulsion between the d-electrons of the impurity. In order to make the problem 
analytically tractable we consider the case of a one impurity problem. 

This problem has been treated using the flow-equation method by Crisan et.al. 
[^for U = and a van-Hove density of states. In the next section we consider the case 
U ^ which is a more realistic case. 
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3 The Flow Equations 



The flow equations method, which is in fact a renormalization procedure apphed in the 
Hamiltonian formalism has been applied in the solid-state theory by the Wegner|jl|, 
Kehrein and Mielke Q and has the main point the diagonalization of the Hamiltonian 
which describes the system by a continuous unitary transformation ry(/) which lead 
to a Hamiltonian. H{1) with the parameters functions of the flow parameter /. This 
transformation satisfies: 

dH{l) 



dl 

where r]{l) can be calculated from: 



[V{1),H{1)] (4) 



rj{l) = [Hoil),Hintil)] (5) 

In order to solve the differential equations using the Hamiltonian. (1) we introduce, 
following 0] the initial values: 



and the renormalized value for I — > oo 



4(0 


= ek{l = 0) 


4(0 


= edil = 0) 




= V{1 = 0) 


U\l) 


= U{1 = 0) 


— > oo 






= efc(oo) 


ef(0 


= ed(oo) 


V^il) 


= V{oo) 


U^{1) 


= U{oo) 



(6) 



(7) 



Using the Eq. (5) and the general method (see Ref.[U|) we calculate r]{l) as: 

r,(/) = r/W (0 + 7?(i) (0 + r/(2) (/) + r]^^) (/) + (;) (g) 

and obtain: 

k,cr 

+ Y.^udtA,. + 4,.d-^) (9) 

k,a- 
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where: 

r)k = i^k - ^d)Vk,d (10) 

and: 

^k,a = -AF_fe(7 (11) 

In Eq.(ll) a =t,i correspond to cr = ±1 in the right side. The higher order con- 
tributions will be given as functions of expressions given by Eqs.(lO) and (11) and 
by: 

&k,a = V-k^O- + ^k,aed (12) 



as: 



Vkl = i^kVk - €dVk - Ae_fe,<,(7) (13) 

= ^k{Vk,aVki + Vki,aVk) - A{^^k,aVk, - ^ki,aV-k)(T (14) 

= {-AV.ka - efcGfc,, - e^e^,,) (15) 

= [HVkraVk + r]kaVk,)a + ek{Ckl,aV-k + C-k,aVk,)] (16) 



,(3) 



{ekVk,aU - edr]k,aU + AC/ C-k,a<^) (17) 
r/^i = {Ari^k,aUa - eu^kU - ea^kU) (18) 
If we take the spin orientation as cr = 1 (the non-magnetic states) the flow equations 



are: 



m ^ (l)r _ , , 

Vk i^k ea + j^^^ _ ^ ^ ^^jj^^ -€, + [/]- 

^ = -4E4^V. (19) 
fc 

dA 1 (2)x. 



iv(o)^''^ '-^ 

where is the Fermi function. 
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4 Solutions of the flow equations 



Using the spectral function 

J{e,l) = Y;^ViS{e-ek{l)) (20) 
k 

and the factorization = Vfc/(efc,0 tfic Eqs. (17) becomes as follows: 

ded f dJ{e,l) [ed-e + Ui][ed + e + U2]-A^ 

dl J ^ dl [ed-e + U][ed + e + U2][ed-e]-A''{ed-e-U) ^ ' 

where: 

Ui = U{l + n{e)) 

U2 = C/(l-n(e)) (22) 

The equation for U becomes: 

dU .f,, dJ{e,l) u[u, + ed + e] 

dl J dl [ed-e + U][ed-e][U2 + ed + e]-A^[ed-e-U] ^ ' 

These equations contain so we have to transform the equation for A as: 

dl N{0)J dl [ed-e + U][ed-e][U2 + ed + e\-A'^[ed-e-U] ^ ' 
In the Eqs. (22),(23),(24)we have ^"^q/^ which is obtained from Eqs. (18) as: 

^ = 2J(M)/(M)l..-. + p;37Tt7Ii^^i;TlFA^> (25) 
and because we have this relation, a supplementary equation for Vk{l) gives no more 
information about the system. The solutions of these equations will be obtained at 
r = (n(e) = 1 — 6(e)) and if we take a concrete form for the J(e, / = 0) as 

21/2 r 

■'(<^.' = 0) = ^ = - P6) 

r = where D is the bandwidth in the limit U » D » we obtained from 
Eqs. (19) using conditions (6), (7): 

e<i = ^d -- — ji 1 - 7T + arctanh^ + ln^\ 

A| = A|(l + /i) (27) 

rrii^ 2r Cd + D ed + U^-D, 
= {7"H [In — + ln- 
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where: 



7, = ' [2ln d +^ +^ +ln i 

^(ef + C^ ^ (28) 



5 Conclusions 

Using the flow equations method we showed that for a BCS superconductor with mag- 
netic correlations describe by the Anderson Hamiltonian with U ^ we calculated the 
energy of impurity the order parameter A and th e energy U. For a large D we get 

Ar = Ai U^ = (29) 

and th energy of the impurity presents a small variation as function of D. 
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